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Abstract 

Some years ago we proposed a new approach for the analysis of galaxy 
and cluster correlations based on the concepts and methods of modern Sta- 
tistical Physics. This led to the surprising result that galaxy correlations 
are fractal and not homogeneous up to the limits of the available cata- 
logues. In the meantime many more redshifts have been measured and 
we have extended our methods also to the analysis of number counts and 
angular catalogues. This has led to a complete analysis of all the available 
data that we are going to present in detail in this lecture. In particular 
we will discuss the properties of the following catalogues: CfA, Perseus- 
Pisces, SSRS, IRAS, Stromlo-APM, LEDA, Las Campanas and ESP for 
galaxies and Abell and ACO for clusters. The result is that galaxy struc- 
tures are highly irregular and self-similar. The usual statistical methods, 
based on the assumption of homogeneity, are therefore inconsistent for 
all the length scales probed until now. A new, more general, conceptual 
framework is necessary to identify the real physical properties of these 
structures. In the range of self-similarity theories should shift from "am- 
plitudes" to "exponents". The new analysis shows that all the available 
data are consistent with each other and show fractal correlations (with di- 
mension D ~ 2) up to the deepest scales probed until now (lOOO/i -1 Mpc) 
and even more as indicated from the new interpretation of the number 
counts. The distribution of visible matter in the universe is therefore 
fractal and not homogeneous. The evidence for this being very strong up 
to 150/1 -1 Mpc due to the statistical robustness of the data and progres- 
sively weaker (statistically) at larger distances due to the limited data. 
In addition the luminosity distribution is correlated with the space distri- 
bution in a specific way characterized by multifractal properties. These 
facts lead to fascinating conceptual implications about our knowledge of 
the universe and to a new scenario for the theoretical challenge in this 
field. 



1 



1 Introduction 



This lecture reports the evidences for the fractal properties of visible matter in 
the universe in relation to the debate with Prof. M. Davis. We have also included 
new elements stimulated by the debate, like the analysis of the Stromlo-APM 
galaxy catalogue and a discussion of the angular data. For a comprehensive and 
detailed report the reader may refer to and to a more recent review Q]. 

From the experimental point of view there are four main facts in Cosmol- 
ogy: - The space distribution of galaxies and clusters: the recent availability 
of several three dimensional samples of galaxies and clusters permits the direct 
characterization of their correlation properties. - The cosmic microwave back- 
ground radiation (CMBR), that shows an extraordinary isotropy and an almost 
perfect black body spectrum. - The linearity of the redshift- distance relation, 
usually known as the Hubble law. This law has been established by measur- 
ing independently the redshift and the distance of galaxies. - The abundance 
of light elements in the universe. Each of these four points provides an inde- 
pendent experimental fact. The objective of a theory should be to provide a 
coherent explanation of all these facts together. 

Our work refers mainly to the first point, the space distribution of galax- 
ies and clusters which, however, is closely related to the interpretation of all 
the other points. In particular we claim that the usual methods of analysis 
are intrinsically inconsistent with respect to the properties of these samples. 
The correct statistical analysis of the experimental data, performed with the 
methods of modern Statistical Physics, shows that the distribution of galaxies 
is fractal up to the deepest observed scales [Q |Q. This result has caused a 
strong debate in the field because it is in contrast with the usual assumption of 
large scale homogeneity which is at the basis of most theories. Actually homo- 
geneity represents much more than a working hypothesis for theory, it is often 
considered as a paradigm or principle and for some authors it is conceptually 
absurd even to question it For other authors instead homogeneity is just the 
simplest working hypothesis and the idea that nature might actually be more 
complex is considered as extremely interesting ||. These two points of view 
are not so different after all because, if something considered absurd becomes 
real, then it is indeed very exciting. Given this situation it may be interesting to 
analyze why this question develops such strong feelings. This will help us to dis- 
tinguish opinions from bare facts and to place the discussion in the appropriate 
perspective. 

Most of theoretical Physics is based on analytical functions and differential 
equations. This implies that structures should be essentially smooth and irreg- 
ularities are treated as single fluctuations or isolated singularities. The study of 
critical phenomena and the development of the Rcnormalization Group (RG) 
theory in the seventies was therefore a major breakthrough |?j |8). One could 
observe and describe phenomena in which intrinsic self-similar irregularities de- 
velop at all scales and fluctuations cannot be described in terms of analytical 



2 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.: 

1.0 
0.8 
0.6 
0.4 
0.2 
0.0 

1.0 i — , — fl n . 

0.8 -C^'V* i% S;:J 
0.6 - , 

0.4 - 1 5 : 
0.2 
0.0 



1.0 



- 10" 



0.0 0.2 0.4 0.6 0.: 

x 



1.0 0.0 0.2 0.4 0.6 0.! 

x 



1.0 



Figure 1: Example of analytical and nonanalytic structures. Top panels (Left) 
A cluster in a homogenous distribution. (Right) Density profile. In this case 
the fluctuation corresponds to an enhancement of a factor 3 with respect to the 
average density. Bottom panels (Left) Fractal distribution in the two dimen- 
sional Euclidean space. (Right) Density profile. In this case the fluctuations 
are non-analytical and there is no reference value, i.e. the average density. The 
average density scales as a power law from any occupied point of the structure. 
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functions. The theoretical methods to describe this situation could not be based 
on ordinary differential equations because self-similarity implies the absence of 
analyticity and the familiar mathematical Physics becomes inapplicable. In 
some sense the RG corresponds to the search of a space in which the problem 
becomes again analytical. This is the space of scale transformations but not the 
real space in which fluctuations are extremely irregular. For a while this peculiar 
situation seemed to be restricted to the specific critical point corresponding to 
the competition between order and disorder. In the past years instead, the de- 
velopment of Fractal Geometry , has allowed us to realize that a large variety 
of structures in nature are intrinsically irregular and self-similar (Fig.|l|). 

Mathematically this situation corresponds to the fact that these structures 
are singular in every point. This property can be now characterized in a quanti- 
tative mathematical way by the fractal dimension and other suitable concepts. 
However, given these subtle properties, it is clear that making a theory for the 
physical origin of these structures is going to be a rather challenging task. This 
is actually the objective of the present activity in the field jl3]. The main dif- 
ference between the popular fractals like coastlines, mountains, trees, clouds, 
lightnings etc. and the self-similarity of critical phenomena is that criticality 
at phase transitions occurs only with an extremely accurate fine tuning of the 
critical parameters involved. In the more familiar structures observed in nature 
instead the fractal properties are self-organized, they develop spontaneously 
from the dynamical process. It is probably in view of this important difference 
that the two fields of critical phenomena and Fractal Geometry have proceeded 
somewhat independently, at least at the beginning. 

The fact that we are traditionally accustomed to think in terms of analytical 
structures has crucial consequences on the type of questions we ask and on the 
methods we use to answer them. If one has never been exposed to the subtleness 
on nonanalytic structures, it is natural that analyticity is not even questioned. 
It is only after the above developments that we could realize that the property of 
analyticity can be tested experimentally and that it may or may not be present 
in a given physical system. 

We can now appreciate how this discussion is directly relevant to Cosmol- 
ogy by considering the question of the Cosmological Principle (CP). It is quite 
reasonable to assume that we are not in a very special point of the universe and 
to consider this as a principle, the CP. The usual mathematical implication of 
this principle is that the universe must be homogeneous. This reasoning implies 
the hidden assumption of analyticity that often is not even mentioned. In fact 
the above reasonable requirement only leads to local isotropy. For an analytical 
structure this also implies homogeneity ||. However, if the structure is not 
analytical, the above reasoning does not hold. For example a fractal structure 
has local isotropy but not homogeneity. In simple terms one observes the same 
mix of structures and voids in different directions (statistical isotropy). This 
means that a fractal structure satisfies the CP in the sense that all the points 
are essentially equivalent (no center or special points) but this does not imply 
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that these points are distributed uniformly. In this respect the saturation of 
the dipole moment cannot be considered as a test of homogeneity because this 
property is also present in fractal structures [pj. 

This important distinction between isotropy and homogeneity has not been 
adequately appreciated in M. Davis' lecture, and it has other important con- 
sequences. For example it clarifies that drawing conclusions about real galaxy 
correlations from the angular distributions alone can be rather misleading as 
we are going to see in detail later. In addition, from this new perspective, the 
isotropy of the CMBR may appear less problematic in relation to the highly 
irregular three dimensional distribution of matter and this may lead to theoret- 
ical approaches of novel type for this problem. In the present work, however, 
we are going to limit our discussion to the optimal way to analyze the data 
provided by the galaxy catalogues from the broader perspective in which ana- 
lyticity and homogeneity are not assumed a priori but they are explicitly tested. 
The main result is that the data become actually consistent with each other and 
coherently point to the same conclusion of fractal correlations up to the present 
observational limits. In this lecture we present a colloquial discussion of the 
main results and their implications. 

2 Statistical Methods and Correlation Proper- 
ties 

- Usual arguments: Before the extensive redshift measurements of the 80s the 
information about the galaxy distribution was only in terms of the two angular 
coordinates. These angular distributions appeared rather smooth at relatively 
large angular scale, like for example the lower part of Fig.|^. Assuming that 
this smoothness corresponds to a real homogenization in 3-d space and esti- 
mating the characteristic depth of the angular catalogue from the magnitudes 
"a characteristic length" r — bh^Mpc has been estimated [B. The idea was 
that beyond such a distance the 3-d galaxy distribution would become gradually 
homogeneous and it could be well approximated by a constant galaxy density. 
This average density, apart from the eventual Dark Matter, would be the one to 
use into Einstein equations to derive the Friedmann metric and the other usual 
concepts. 

Later, the measurements of the galaxy redshifts plus the Hubble law provided 
also the absolute distances and could identify the position of galaxies in space. 
However the 3-d galaxy distributions turned out to be much more irregular with 
respect to their angular projections and showed large structures and large voids, 
as shown in the upper part of Figj^. At first these irregular structures appeared 
to be in contradiction with the picture derived from the angular catalogues and, 
as we are going to see, they really are. However in 1983 a correlation analysis of 
the 3-d distribution (CfAl catalog |Q) was performed by Davis and Peebles [jl3| 
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Figure 2: A slice of the three dimensional galaxy distribution (a) (b) compared 
with the corresponding (c) angular distribution (the portion between 26.5° and 
32.5° - [11]). Note that the angular distribution appears relatively homogenous 
while the real three distribution in space is much more irregular. In particu- 
lar this picture shows the so-called Great Wall which extends over the entire 
sample (at least 170h~ 1 Mpc). Note that it is impossible to identify structures 
of this typo from the angular data alone. We also show the size of the galaxy 
" correlation- length" (r = 5h~ 1 Mpc) derived from the standard analysis. The 
more general analysis we discuss here implies instead that an eventual correla- 
tion length should be larger than any observable structure, i.e. 3> 170h~ 1 Mpc 
and that the present data show well defined fractal properties up to the sample 
limits. These pictures clarify the intrinsic ambiguity of the angular catalogs in 
relation to the arguments of Prof. Davis. They show that a fractal distribution 
in 3-d (top) projects in a smooth angular distribution (bottom). No one could 
have predicted the real 3-d properties from the angular distribution alone. 
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and the result was again that the correlation length was ro = 5h~ 1 Mpc as for the 
angular catalogues. This seemed to resolve the puzzle because it was interpreted 
as if a relatively small correlation length can be consistent with the observation 
of large structures. This value for ro has not been seriously questioned even 
after the observation of huge structures, like the galaxy wall, that extend up 
to - Y7Qh- 1 Mpc. Fig| is a clear example of the smoothing effect of angular 
projections, and it already gives a clear indication of the compatibility of a 
fractal structure in 3-d with a smoother projection. 

The usual correlation analysis is performed by estimating at which distance 
(ro) the density fluctuations are comparable to the average density in the sample 
(£(r) = (n(O)n(r)) / (n) 2 — 1; £(r ) = 1). Now everybody agrees that there are 
fractal correlations at least at small scales. The important physical question is 
therefore to identify the distance Ao at which, possibly, the fractal distribution 
has a crossover into a homogeneous one. This would be the real correlation 
length beyond which the distribution can be approximated by an average den- 
sity. The problem is therefore to understand the relation between ro and Ao: are 
they the same or closely related or do they correspond to different properties? 
This is actually a subtle point with respect to the concepts discussed in the 
introduction. In fact, if the galaxy distribution becomes really homogeneous at 
a scale Ao within the sample in question, then the value of ro is related to the 
real correlation properties of the distribution and one has ro ~ Ao/2. If, on the 
other hand, the fractal correlations extend up to the sample limits, then the 
resulting value of ro has nothing to do with the real properties of the galaxy 
distribution but it is fixed just by the size of the sample |Q. 

New Perspective: Given this situation of ambiguity with respect to the real 
meaning of ro it is clear that the usual correlation study in terms of the func- 
tion £(r) is not the appropriate method to clarify these basic questions. The 
essential problem is that, by using the function £(r), one defines the amplitude 
of the density fluctuations by normalizing them to the average density of the 
sample in question. This implies that the observed density should be the real 
one and it should not depend on the given sample or on its size apart from 
Poisson fluctuations. However, if the distribution shows long range (fractal) 
correlations, this approach becomes meaningless. For example if one studies a 
fractal distribution with £(r) a characteristic length ro will be identified, but 
this is clearly an artifact because the structure is characterized exactly by the 
absence of any defined length ffl. 

The appropriate analysis of pair correlations should therefore be performed 
using methods that can check homogeneity or fractal properties without assum- 
ing a priori either one. The simplest method to do this is to consider directly 
the conditional density T(r) ~ (n(O)n(r)} without comparing it to the aver- 
age. There are several other methods that have been discussed elsewhere [Q 
[||. This is not all however because one has also to be careful not to make 
hidden assumptions of homogeneity in the specific procedure to evaluate these 
correlations. For example, if a galaxy is close to the boundary of the sample, it 
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is possible that the sphere of radius r around it, where the conditional density is 
computed, may lie in part outside the sample boundary. In this case the usual 
procedure is to use weighting schemes of various types to include also these 
points in the statistics. In this way one implicitly assumes that the fraction of 
sphere contained in the sample is sufficient to estimate the properties of the full 
sphere. This implies that the properties of a small volume are assumed to be 
the same as for a larger volume (the full sphere). This is a hidden assumption 
of homogeneity that should be avoided by including only the properties of those 
points for which a surrounding sphere of radius r is fully included in the sam- 
ple. These procedures are fully standard in modern statistical mechanics and a 
detailed description can be found in §. This means that the statistical validity 
of a sample is limited to the radius of the largest sphere that can be contained 
in the sample. We call this distance R s and it should not be confused with the 
sample depth Rd, which can be in general much larger, depending on the survey 
geometry. 

In 1988 we reanalyzed the CfAl catalogue jl5|. The result was that the 
catalogue has statistical validity up to R s — 20h~ 1 Mpc and, up this length, 
it shows well defined fractal correlations. This shows therefore that the "cor- 
relation length" ro = 5h~ 1 Mpc derived by |l3| was a spurious result due to 
an inappropriate method of analysis and it has nothing to do with the real 
correlation properties of the system. A similar analysis of the Abell cluster 
catalogue also showed fractal properties up to R s — 80h~ 1 Mpc so that also 
the cluster "correlation length" rg = 25h~ 1 Mpc |ll| should be considered as 
spurious. One consequence of these results was that the so called galaxy-cluster 
mismatch could be automatically eliminated by the appropriate analysis. Also 
other properties like 6N/N, directly related to ro, suffer from the same con- 
sistency problems because of the lack of a reference value [|l7| . This situation 
led to a rather controversial debate in the field. In the meantime many more 
data have became available and we have performed a complete analysis of all 
the data for galaxies and clusters. In the following we report the main results. 

3 Analysis of the Galaxy Distributions 

Here we discuss the correlation properties of the galaxy distributions in terms of 
volume limited catalogues jl|] arising from most of the 50.000 redshift measure- 
ments that have been made to date. A first important result will be that the 
samples are statistically rather good and their properties are in agreement with 
each other. This gives a new perspective because, using the standard methods 
of analysis, the properties of different samples appear contradictory with each 
other and often this is considered to be a problem of the data (unfair samples) 
while, we show that this is due to the inappropriate methods of analysis. In 
addition essentially all the catalogues show well defined fractal correlations up 
to their limits and the fractal dimension is D ~ 2. The few exceptions to this 
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Table 1: The volume limited catalogues are characterized by the following pa- 
rameters: - Rdihr 1 Mpc) is the depth of the catalogue - is the solid angle - 
-R s (/i _1 Mpc) is the radius of the largest sphere that can be contained in the 
catalogue volume. This gives the limit of statistical validity of the sample. - 
r {h~ 1 Mpc) is the length at which £(r) = 1. - A is the eventual real crossover 
to a homogeneous distribution that is actually never observed. The value of rg is 
the one obtained in the deepest sample. The CfA2 and SSRS2 data are not yet 
available. See Ref [13], [17]-[23] and [31] (distance are expressed in h~ 1 Mpc). 



Sample 


n (sr) 


Rd 


R s 


J*0 


D 


A 


CfAl 


1.83 


80 


20 


6 


1.7 ±0.2 


> 80 


CfA2 


1.23 


130 


30 


10 


2.0 


? 


PP 


0.9 


130 


30 


10 


2.0 ±0.1 


> 130 


SSRS1 


1.75 


120 


35 


12 


2.0 ±0.1 


> 120 


SSRS2 


1.13 


150 


50 


15 


2.0 


? 


Stromlo-APM 


1.3 


100 


30 


10 


2.2 ±0.1 


> 150 


LEDA 


47T 


300 


150 


45 


2.1 ±0.2 


> 150 


LCRS 


0.12 


500 


18 


6 


1.8 ±0.2 


> 500 


IRAS 1.2 J y 


47T 


80 


40 


4.5 


2.0 ±0.1 


~ 25 


ESP 


0.006 


700 


10 


5 


1.9 ±0.2 


> 800 



result will be discussed and interpreted in detail. 

The main data of our correlation analysis are collected in Fig. | (left part) 
in which we report the conditional density as a function of scale for the vari- 
ous catalogues. The relative position of the various lines is not arbitrary but 
it is fixed by the luminosity function, a part for the cases of IRAS and SSRS1 
for which this is not possible. The properties derived from different catalogues 
are compatible with each other and show a power law decay for the conditional 
density from lh Mpc to 150h Mpc without any tendency towards homog- 
enization (flattening). This implies necessarily that the value of tq (derived 
from the £(r) approach) will scale with the sample size R s as shown also from 
the specific data about tq of the various catalogues jl5| p2| . The behaviour 
observed corresponds to a fractal structure with dimension D ~ 2. The smaller 
value of CfAl was due to its limited size. An homogeneous distribution would 
correspond to a flattening of the conditional density which is never observed. 
Usually the £(r) approach also leads to a smaller value of D (or a larger value 
of 7 = 3 — D) as derived from the small scale properties. This is due to the fact 
that the fit is made close to ro and it is affected by the fact that £(r) (in log 
coordinates) is becoming steeper because it is approaching the value zero. 

It is interesting to compare the analysis of Fig.0 with the usual one, made 
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Figure 3: Full correlation for the various available redshift catalogues in the 
range of distances 0.1 -j- 1000h~ 1 Mpc. A reference line with a slope —1 is 
also shown (i.e. fractal dimension D = 2). Up to ~ 150h~ 1 Mpc the density 
is computed by the full correlation analysis, while above ~ 150h~ 1 Mpc it is 
computed through the radial density. For the full correlation the data of the 
various catalogues are normalized with the luminosity function and they match 
very well with each other [2] . This is an important test of the statistical validity 
and consistency of the various data. In the insert panel it is shown the schematic 
behavior of the radial density versus distance computed from the vertex (see 
text) [27]. The behaviour of the radial density allows us to extend the power 
law correlation up to ~ 1000/i _1 Mpc. However a rescaling is necessary to match 
the radial density to the conditional density [2] . 
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Figure 4: Usual analysis based on the function £(r) of the same galaxy cata- 
logues of Fig. 3. This analysis is based on the a priori and untested assumption 
of the analyticity and homogeneity. These properties are not present in the real 
galaxy distributions and the results appear therefore rather confusing. This lead 
to the impression that galaxy catalogues are not good enough and to a variety of 
theoretical problems like the galaxy cluster mismatch, luminosity segregation, 
the linear and non linear evolution, etc.. The situation changes completely and 
it becomes rather clear if one adopts the more general framework that is at the 
basis of Fig. 3. In the insert panel we show the dependence of r on R s for all the 
catalogs. The linear behaviour is a consequence of the fractal nature of galaxy 
distribution in these samples [2]. 
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with the function £(r), for the same galaxy catalogues. This is reported in Fig.[| 
and, from this point of view, the various data appear to be in strong disagree- 
ment with the each other. This is due to the fact that the usual analysis looks 
at the data from the perspective of analyticity and large scale homogeneity 
(within each sample). These properties are never tested and they are actually 
not present in the real galaxy distributions so the result is rather confusing 
(Fig.|). Once the same data are analyzed within a broader perspective the 
situation becomes clear (Fig.||) and the data of different catalogues result in 
agreement with each other. In addition in the insert of Fig.^j we show the de- 
pendence of ro on R s for all the catalogs. The linear behaviour is a consequence 
of the correlation properties of Fig.|| and it provides an additional evidence of 
fractal behaviour to all scales. In this respect, the proposed luminosity bias ef- 
fect mentioned by M. Davis appears essentially irrelevant while, on the contrary, 
the linear dependence of ro on R s is very clear. 

The new picture allows us to make clear predictions for the value of ro for the 
forthcoming catalogs CfA2 and SLOAN. Considering the predicted behaviour of 
ro ~ 0.3i? s (for D ~ 2 - see insert of Fig.|]) we expect that in CfA2 one should 
have r w 15~20h~ 1 Mpc while for the full SLOAN catalog r » m+Wh^Mpc 
(of course these values depend on the solid angle of the survey and the volume 
limited sample considered in the way precisely discussed previously). We stress 
however that these predictions refer to the full solid angle catalogs and for 
subsamples one should consider the corresponding value of R s . 

It is important to remark that analyses like those of Fig.|| have had a pro- 
found influence in the field in various ways: first the various catalogues appear 
in conflict with each other. This has generated the concept of " not fair sample" 
and a strong mutual criticism about the validity of the data between different 
groups. In other cases the discrepancy observed in Fig.^| have been considered 
as real physical problems for which various theoretical approaches have been 
proposed. These problems are, for example, the galaxy-cluster mismatch, lumi- 
nosity segregation, the richness clustering relation and the linear and non linear 
evolution of the perturbations corresponding to the " smaW or " large" ampli- 
tudes of fluctuations. We can now see that all this problematic is not real and it 
arises only from a statistical analysis based on inappropriate and to restrictive 
assumptions that do not find a correspondence in the physical reality. It is also 
important to note that, even if the galaxy distribution would eventually become 
homogeneous at some large scale, the use of the above statistical concepts is 
anyhow inappropriate for the range of scales in which the system shows fractal 
correlations as those shown in Fig.[| 

Contrary to the claims of Prof. Davis we would like to stress that a fractal 
distribution has a very strong property: it shows global power-law correlations 
up to the sample depth. Such correlations cannot be due neither to an inhomo- 
geneous sampling of an homogeneous distribution, nor to some selection effects 
that may occur in the observations. Namely, suppose that a certain kind of 
sampling reduces the number of galaxies as a function of distance. Such an ef- 
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feet, in no way can lead to long range correlations, because when one computes 
r(r), one makes an average over all the points inside the survey. In any case this 
possible bias could be detected by a difference in the values of T(r) at different 
depths, which is not observed 0. We observe instead that all the catalogues, 
independently on their completeness, show precisely the same correlation prop- 
erties. 

In connection with the conditional density decay of Fig|], it is interesting 
to note that the Hubble law (redshift versus distance) has been experimentally 
tested in the same range of scales. These two experimental facts show that 
the Hubble law is compatible with a fractal universe, contrary to the usual 
theoretical interpretation |l7[] p^| . 

In FigJ^ we report the power spectrum analysis (Fourier conjugate of of 
some catalogues. Also in this for £(r), there is a specific dependence of 

the properties on the sample size, in full agreement with the direct correlation 
analysis of Fig.|[ In particular, for a fractal structure, the amplitude of the 
power spectrum is a function of the sample size and the shape is characterized 
by a turnover: both these features, bending and scaling, are a manifestation of 
the finiteness of the survey volume, and cannot be interpreted as the convergence 
to homogeneity, nor to a power spectrum flattening. A detailed discussion of 
the power spectrum can be found in [ p7| . 

Essentially similar results have been obtained for the Abell and ACO cluster 
catalogues § g @. Also in these cases we observe a power law behaviour 
for the cluster correlations with D ~ 2 and without any tendency towards 
homogenization up to ~ 80h~ 1 Mpc. 

All these results imply that the previous "correlation lengths" of 5h~ 1 Mpc 
and 25h~ 1 Mpc, introduced for galaxies and clusters, are spurious and no real 
correlation length can be defined from the data. Therefore the much discussed 
mismatch between galaxy and cluster correlations, that is also at the basis 
of various theories for structure formation, does not actually exist. Cluster 
correlations correspond just to the continuation of galaxy correlations at larger 
scales. In the language of Statistical Physics, cluster catalogues are the coarse 
grained version of galaxy catalogues. 

4 Statistical Validity of Catalogues 

Often the concept of fair sample has been used to mean homogeneous sample. 
For this reason various samples are declared as not fair just because they are not 
homogeneous. We have seen instead that it was the method of analysis to be 
inappropriate while the samples are actually rather good in a statistical sense. 
This means that one can derive from them a statistically valid information about 
their correlation properties. 

In relation to the statistical validity it is interesting to consider the IRAS 
catalogues because they seem to differ from all the other ones and to show some 



13 




a cede fMpc/li') 



Figure 5: Power spectrum (PS) of several samples of LEDA together with the 
CfA2 power spectra (from [22] The amplitude of the PS depends on the sample 
depth, and it is larger in the case of LEDA, because in this case the effective 
depth is larger than the CfA2 one (see Table 1). The curved shape of the PS 
is also spurious that does not correspond to any intrinsic feature in the case of 
fractals, but it is only due to the finiteness of the sample volume [25]. We also 
plot a CDM-like fit both to the CFA2 data and the LEDA data. On the right, 
we plot the COBE data. No reasonable CDM curve can fit the data and COBE 
(from [25]) 
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tendency towards homogenization at a relatively small scale Actually the point 
of apparent homogeneity is only present in some samples, it varies from sample 
to sample between ~ 15 -j- 25h~ 1 Mpc and it is strongly dependent on the dilu- 
tion of the sample. Considering that structures and voids are much larger than 
this scale and that the IRAS galaxies appear to be just where luminous galaxies 
are it is clear that this tendency appears suspicious. One of the characteristic 
of the IRAS catalogues with respect to all the other ones is an extreme degree 
of dilution: this catalogue contains only a very small fraction of all the galaxies. 
It is important therefore to study what happens to the properties of a given 
sample if one dilutes randomly the galaxy distribution up to the IRAS limits. 
A good test can be done by considering the Perseus Pisces catalogue and elim- 
inating galaxies from it. The original distribution shows a well defined fractal 
behavior. By diluting it to the level of IRAS one observes an artificial flattening 
of the correlations 28 (see Fig.||). This effect does not correspond to a real 
homogenization but it is due to the of dilution. In fact it can be shown that 
when the dilution is such that the average distance between galaxies becomes 
comparable with the largest voids (lacunarity) of the original structure there is 
a loss of correlation and the shot noise of the sparse sampling overcomes the 
real correlations and produces an apparent trend to homogenization [27]. This 
allows us to reconcile this peculiarity of the IRAS data with the properties of 
all the other catalogues. Analogous considerations for other sparse samples like 
QDOT and the Stromlo-APM samples @]. 

The correlations discussed up to now are well defined statistically but limited 
to the radius of the largest sphere that can be contained in the sample R s . For 
example for Las Campanas the depth is very large Rd ~ 500h~ 1 Mpc but R s is 
only 20h~ 1 Mpc because the sample is very thin. So, it is not surprising that 
the value of ro is also small (6 -r- 7h~ x Mpc). Given this situation it would be 
very interesting to find some method that is limited by Rd instead of R s . 

Galaxy samples have typically a conic shape and, in this respect, they are 
three dimensional objects. Considering a volume limited sample, if one simply 
counts the number of galaxies within a distance r from the vertex this number 
should go like r 3 for a homogeneous distribution and like r D for a fractal one. In 
order to compare with the previous correlation analysis it is actually convenient 
to consider the density instead of the total number. The problem with such an 
analysis is that one cannot average from different observational points but the 
count is from a single point, the vertex. This situation corresponds to a reduced 
statistical quality that should be carefully analyzed. 

At very small distances we are not going to find any point. When a distance 
of the order of the minimal one between galaxies is reached, we begin to have a 
signal but this is strongly affected by finite size effects as shown schematically 
in the insert of FigJ^. Eventually at some distance A, the number of points 
becomes large enough that one reaches the correct scaling behavior. It can 
be shown that this characteristic length for the statistical significance of this 
method is about ten times the minimal distance between galaxies. For various 
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Figure 6: Top panels: The conditional density (left) for a volume limited 
sample of the full Perseus-Pisces redshift survey (right). The percentage of 
galaxies present in the sample is ~ 6%. The slope is —7 = — 1. Bottom panels: 
In this case the percentage of galaxies is ~ 1%, and the number of galaxies 
is the same of the IRAS 1.2 Jy sample in the same region of the sky. We can 
see that at small scale we have a 1/r 3 decay just due to the sparseness of the 
sample, while at large scale the shot noise of the sparse sampling overcomes the 
real correlations and produces an apparent trend to homogenization [27] . In our 
opinion, this effect, due to sparsness of this sample, is the origin of the apparent 
trend towards homogencization observed in some of the IRAS samples. 
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catalogues the value of A is appreciably smaller that Rd and, in these cases, 
a useful information can be obtained for the length scales between A and Rd- 
A detailed discussion of this method can be found in p8| |. This approach is 
quite useful because it allows us to use the thin deep catalogues up to their 
total depth Rd- Particularly interesting in this respect is the ESP catalogue 
whose depth extends up to 900h~ 1 Mpc. Also for the Las Campanas sample 
it is possible to obtain some information despite the peculiar and unfortunate 
luminosity selections of this catalogue. The results of these deep catalogue are 
reported in Fig.|| together with those of the other catalogues discussed before. 
The amplitude of the density computed from the vertex for all the catalogues 
is systematically shifted by about a factor of 3 (that has been rescaled in Fig.||) 
with respect to the full correlation analysis. This shift is due to the fact that 
usually, observations point towards zones that are rich of galaxies. A detailed 
discussion of this effect is reported in Q . 

The behaviour of the density decay from the vertex shows the same power 
law behavior of the full correlation analysis of Fig.|| (left part) but this property 
is now shown to extend up to 900h~ 1 Mpc. The agreement between different 
catalogues and different methods of analysis is remarkable and it shows a well 
defined fractal behavior for the galaxy correlations extending from 1 to about 
W00h~ 1 Mpc (Fig|j). This analysis refutes therefore the early comments about 
a possible homogenization in Las Campanas, based on the visual impression of 
the data. 

In relation to the lecture by Prof. Davis, it should be noted that the function 
n(z) (redshift counts in a magnitude limited sample) it can be shown Q that 
the behaviour of such a quantity is mostly related to the luminosity selection 
function of the survey rather than to the behaviour of the space density. In 
particular, such a quantity cannot be used to distinguish between fractal or 
homogeneous properties of the galaxy distribution. 

5 Number Counts and Angular Correlations 

The above discussion of the density decay from the vertex brings us naturally 
to the problem of the galaxy number counts that is also performed in this 
way, namely by counting from the origin. There are however some relevant 
differences because all the properties discussed up to here refer to volume limited 
samples while the galaxy counts are defined by their apparent magnitude. The 
behaviour of the number versus magnitude relation (N(< m)) is reported in 
Fig.^. The exponent a defined by this relation can be related to the dimension 
D of the galaxy distribution (a — D/h). For small magnitudes (small scales), 
one observes a ~ 0.6 ±0.1 that is usually interpreted as evidence of homogeneity 
(D « 2). At larger scales the value of a decreases to about 0.4 (corresponding 
to D ~ 2) and this deviation is considered as due to galaxy evolution and solid 
angle effects due to the expansion. However, from the direct correlation analysis 
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Figure 7: The galaxy number counts in the B-band, from several surveys. In 
the range 12 < m < 19 the counts show an exponent a ~ 0.6 ± 0.1, while in 
the range 19 < m < 28 the exponent is a ~ 0.4. The amplitude of the galaxy 
number counts for m > 19 (solid line) is computed from the determination of 
the prefactor B of the density n(r) — Bt d ~ 3 (with D = 2 - see text) at small 
scale and from the knowledge of the galaxy luminosity function. The distance is 
computed for a galaxy with M = — 16 and we have used H Q = 7hkmsec~ 1 Mpc 
[27] [2]. 
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in 3-d space, we know that up to some rather large scale the galaxy distribution 
is certainly fractal. Therefore the usual interpretation of a ~ 0.6 ±0.1 at small 
scales in terms of a homogeneous distribution cannot be correct. In this respect 
the insert of Figj^ allows us to clarify the situation. In fact, at relatively small 
scales, the density first raises and then falls when it reaches the correct scaling 
regime. This behaviour can roughly resemble a constant density, especially for 
an integrated quantity. A variety of tests for various catalogues with fractal 
correlations in 3-d, like the Perseus Pisces survey and others, show that the 
number counts give indeed an exponent a ~ 0.6 ±0.1 if the sample is dominated 
by small scales finite size fluctuations. On the other hand, if one computes the 
number counts by performing an average over all the different observers in a 
redshift sample, the value a rs 0.4 is readily found, in agreement with the space 
correlation analysis. In the deep surveys (m > 19) the small scale fluctuations 
are self-averaging so that one can recover the correct properties (i.e. a ~ 0.4) 
without performing any average. 

This implies a completely different interpretation of the number counts. At 
small scales the value a w 0.6 is due to finite size effects and not to a real homo- 
geneity, while at larger scales the value 0.4 corresponds to the correct correlation 
properties of the sample. This implies that galaxy evolution, modification of the 
Euclidean geometry and the K-corrections are not very relevant in the range of 
the present data. In addition the fact that the exponent 0.4 holds up to mag- 
nitudes 27 -j- 28 seems to indicate that the fractal properties may extend up to 
2 -j- 3000h~ 1 Mpc. The amplitude of the galaxy number counts for m > 19 (see 
the solid line Fig.|^) is computed from the determination of the prefactor B of 
the density n(r) = Br D ~ 3 at small scale and from the knowledge of the galaxy 
luminosity function p8| . Quite a remarkable fact if one considers that the Hub- 
ble radius of the universe is supposed to be 4000/i~ 1 Afpc. Such a behaviour 
(a ~ 0.4, i.e. D ~ 2) can be found for galaxies in the different Photometric 
bands, as well as for other astrophysical objects 

We have now all the elements to reinterpret also the angular catalogues. 
These catalogues are qualitatively inferior to the 3-d ones because they only 
correspond to the angular projection and do not contain the third coordinate. 
However the fact that they contain more points with respect to the 3-d cata- 
logues has led some authors to assign an excessive importance to these cata- 
logues (once we have the 3-d ones). In fact, even for a very large number of 
points [N — > oo) the angular distribution is intrinsically degenerate, in the sense 
that it can originate from different 3-d distributions. Actually the interpreta- 
tion of the angular catalogues is quite delicate and ambiguous for a variety of 
reasons: 

- Unlike orthogonal projections, angular projections mix different length 
scales and this gives an artificial randomization of the points. This can be 
illustrated as follows: consider a small solid angle dQ, of an angular catalog. 
This will contain the projection of all the points in the corresponding cone of 
depth Rd- This cone is a three dimensional object, so the intersection with a 
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fractal of dimension D will also be a set of dimension D. Therefore the number 
of points that project in dfl is N(dft, Rd) ~ dSIR® . This implies that for a large 
enough value of Rd any small solid angle dfl will contain some points. This is 
why structures and voids are smeared by the angular projection. This implies 
that the angular projection of a fractal structure can be really homogeneous at 
relatively large angles |2{| Clearly this is an artificial effect and from a 
smooth angular projection one cannot deduce whether the real distribution is 
also smooth. An example of this effect is given by Fig. 2 in which the angular 
projection appears relatively smooth, while the real distribution is much more 
irregular. In fact the large structures and large voids have been identified with 
the redshift catalogues and could not be predicted from the angular data alone. 

- An additional effect is the one due to the finite size effects discussed before 
in case of a single observational point that gives an additional artificial effect of 
homogenization. The so called rescaling of the angular correlations can be un- 
derstood in detail by these effects and it can be shown that the same properties 
can be observed in a fractal distribution pjfl . 

Concerning the debate with M. Davis and the subsequent discussion with J. 
Peebles it should be noted that the angular projections correspond to a complex 
convolution with the magnitude limit and dilution effects in the case of IRAS. 
Contrary to the claim of M. Davis we have already analyzed the propreties of 
angular projections [^9| [Q Q and identified a number of problems, like those 
mentioned here, that were not known before and that make these catalogues 
intrinsically unambiguous. For these reasons we decided to focus on 3-d dis- 
tributions that lead to the consistent and ambiguous results shown in Fig.^. 
Concerning the proposal of Davis and Peebles to generate a fractal distribution 
in 3-d that corresponds to the observed angular projections, this requires the 
tuning of various other properties in addition to the fractal dimension. In par- 
ticular lacunarity, morphology, magnitude limits and dilution effects all play a 
crucial role in the projection and should be tuned to those of real galaxies. In 
any case 3-d galaxy distributions are fractals and their properties do not depend 
on our ability to make this exercise. 



6 Luminosity and Space Distribution 

Up to now we have discussed galaxy correlations only in terms of the set of points 
corresponding to their position in space. Galaxies can be also characterized by 
their luminosity (related to their mass) and the full luminosity distribution is 
then a full distribution and not a simple set. It is natural then to consider the 
possible scale invariant properties of this distribution. This requires a gener- 
alization of the fractal dimension and the use of the concept of multifractality 
[Q. In this language the fractal set discussed until now represents the support 
of the full luminosity distribution. 

A multifractal analysis shows that also the full distribution is scale invariant 
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|§ and this leads to a new and important relation between the Schechter 
luminosity distribution and the space correlation properties. This allows us 
to understand various morphological features (like the fact that large elliptic 
galaxies are typically located in large clusters) in terms of multifractal expo- 
nents. This leads also to a new interpretation of what has been called the 
luminosity segregation effect |$cj. 

The observation that the most luminous elliptical galaxies lie in the core 
of rich galaxy clusters is, in our analysis, a manifestation of the multifractal 
properties, i.e. the self-similar distribution of matter including galaxy masses 
(luminosities). However this observation together with the shift of vq with sam- 
ple depth (interpreted as luminosity limit) has lead various authors Q 
[ [32"| to formulate the qualitative hypothesis that the homogenization crossover 
(related to ro) should be different for galaxies of different luminosity. The fact 
that large voids (> 50h~ 1 Mpc) are empty of galaxies of any type is already a 
disproof of this hypothesis. In addition we have shown that the appropriate 
correlation analysis shows a power law behaviour at any observable scale. This 
implies unavoidably that ro must scale with the sample size because the system 
is self-similar. Most of the analysis of the luminosity segregation effects usually 
do not address the fundamental question whether ro is a physical meaningful 
quantity, that should be addressed with the conditional density. Only once a 
crossover towards homogeneity would be observed in the conditional density 
then "luminosity segregation" questions could eventually be posed in the usual 
terms. 



7 Conclusions and Theoretical Implications 

In summary our main points are: 

• The highly irregular galaxy distributions with large structures and voids 
strongly point to a new statistical approach in which the existence of a 
well defined average density is not assumed a priori and the possibility of 
non analytical properties should be addressed specifically. 

• The new approach for the study of galaxy correlations in all the available 
catalogues shows that their properties are actually compatible with each 
other and they are statistically valid samples. The severe discrepancies 
between different catalogues that have led various authors to consider 
these catalogues as not fair, were due to the inappropriate methods of 
analysis. 

• The correct two point correlation analysis shows well defined fractal corre- 
lations up to the present observational limits, from 1 to 1000/i _1 Mpc with 
fractal dimension D ~ 2. Of course the statistical quality and solidity of 
the results is stronger up to 100^200/i _1 Afpc and weaker for larger scales 
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due to the limited data. It is remarkable, however, that at these larger 
scales one observes exactly the continuation of the correlation properties 
of the small and intermediate scales. 

• These new methods have been extended also to the analysis of the number 
counts and the angular catalogues which are shown to be fully compatible 
with the direct space correlation analysis. The new analysis of the number 
counts suggests that fractal correlations may extend also to scales larger 
that 1000h- 1 Mpc. 

• The inclusion of the galaxy luminosity (mass) leads to a distribution which 
is shown to have well defined multifractal properties. This leads to a 
new, important relation between the luminosity function and that galaxy 
correlations in space. 

• New perspective on old arguments. On the light of these results we can 
now take a standard reference volume in the field (i.e. Peebles 1993 ||) 
and consider the usual arguments invoked for homogeneity from a new 
point of view. These arguments are: (a) number counts: we have seen 
in Sec. 5 that the small scale exponent of number counts is certainly not 
related to homogeneity but to small scale fluctuations. The real exponent 
of the number counts is instead the lower one (i.e. a « 0.4) that indeed 
corresponds to the three dimensional (i.c fractal with D ~ 2) correlation 
properties, (b) SN/N is small for various observations. This point is 
exactly the same as the fact ro is a spurious length. In the absence of a 
reference average one cannot talk about "large" or "small" amplitude of 
fluctuations. In addition, for any distribution, even a fractal one SN/N is 
always small for sizes comparable to the total sample because the average 
is computed from the sample itself, (c) Angular correlations. We have 
seen in Sec. 5 that angular correlations are ambiguous in two respects: first 
the angular projection of a fractal is really uniform at large angles due to 
projection effects, second the angular data are strongly affected by the fi- 
nite size fluctuations that provide an additional artificial homogcnization, 
as in the case of the number counts. The inclusion of these effects recon- 
ciles quite naturally the angular catalogs with the fractal properties in the 
three dimensional ones, (d) X-ray background. The argument that SN/N 
becomes very small for the X-ray background combines the two problems 
discusses before: angular projections and reference average. This angu- 
lar uniformity is analogues, for example, to the Lick angular sample, and 
certainly is not a proof of real homogeneity. 

Finally one should note that all these arguments are indirect and always 
require an interpretation based on some assumptions. The most direct 
evidence for the properties of galaxy distribution arises from the correct 
correlation analysis of the 3-d volume limited samples that has been the 
central point of our work. 
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Theoretical Implications 

From the theoretical point of view the fact that we have a situation char- 
acterized by self-similar structures implies that we should not use concept like 
£(r), ro, 5N/N and certain properties of the power spectrum, because they are 
not suitable to represent the real properties of the observed structures. In this 
respect also the N-body simulations should be considered from a new perspec- 
tive. One cannot talk about "small" or "large" amplitudes for a self-similar 
structure because of the lack of a reference value like the average density. The 
Physics should shift from "amplitudes" towards "exponent" and the methods 
of modern statistical Physics should be adopted. This requires the development 
of constructive interactions between two fields. 

Possible Crossover. We cannot exclude of course, that visible matter may 
really become homogenous at some large scale not yet observed. Even if this 
would happen the best way to identify the eventual crossover is by using the 
methods we have described and not the usual ones. From a theoretical point 
of view the range of fractal fluctuations, extending at least over three decades 
(1 -7- 1000h~ 1 Mpc) 1 should anyhow be addressed with the new theoretical con- 
cepts. Then one should study the (eventual) crossover to homogeneity as an 
additional problem. For the moment, however, no tendency to such a crossover 
is detectable from the experimental data and it may be reasonable to consider 
also more radical theoretical frameworks in which homogenization may simply 
not exist at any scale, at least for luminous matter. 

Dark Matter. All our discussion refers to luminous matter. It would be 
nice if the new picture for the visible universe could reduce, to some extent, 
the importance of Dark matter in the theoretical framework. At the moment 
however this is not clear. We have two possible situations: (i) if Dark matter 
is essentially associated to luminous matter, then the use of FRW metric is not 
justified anymore. This does not necessarily imply that there is no expansion or 
no Big Bang. It implies, however, that these phenomena should be described by 
more complex models, (ii) If Dark matter is homogenous and luminous matter 
is fractal then, at large scale, Dark matter will dominate the gravity field and 
the FRW metric is again valid. The visible matter however remains self-similar 
and non analytical and it still requires the new theoretical methods mentioned 
before. However, this perspective seems to be in contrast with the usual role of 
Dark Matter, that is to generate large potential wells for the initially smoother 
baryonic matter. 

Predictions and Bets 

After the debate with Prof. M. Davis, we tried to assess the predictions 
of the two points of view and we also agreed to make a bet. According to 
the arguments of Prof. M. Davis, the length r — 5h~ 1 Mpc characterizes the 



23 



physical properties of galaxy distributions. Therefore deeper samples like CfA2 
and SLOAN should simply reduce the error bar, which is now about consid- 
ered to be 10%. A possible variation of ro with absolute magnitude, due to a 
luminosity bias (see Sec. 6), is considered plausible but it has never been quan- 
tified. This should be checked by varying independently absolute magnitude 
and depth of the volume limited samples. However, from this interpretation, 
the value of r = hh^Mpc, corresponding to a volume limited of CfAl with 
M = —19.5, should not change when considering in CfA2 and SLOAN volume 
limited samples with the same solid angle Q and the same absolute magnitude 
limit (M = -19.5). 

In our interpretation, instead, ro is spurious, and it scales linearly with the 
radius R s of the largest sphere fully contained in the volume limited samples. 
Therefore we predict for the volume limited sample of CfA2 with M = —19.5 
(with a solid angle of fi = 1.1 sr |23 ) ro « 7h Mpc (if, in the final version 
of the survey the solid angle will be fl = 1.8, the value of R s will increase 
accordingly, and the value of ro will shift up to ~ 9h~ 1 Mpc). Note however 
that for the deepest volume limited CfA2 sample (M < —20) we predict instead 
ro ~ 15 -7- 20h~ 1 Mpc. For the volume limited sample of the full SLOAN with 
M = —19.5 (f2 = 7r), our prediction is that ro ~ 65h~ 1 Mpc. It is clear that 
however, the first SLOAN slice will give smaller values because the solid angle 
will be small. 

About the respective predictions for the full SLOAN (O = tt sr) one of us 
(L.P.) made a bet with Prof. M. Davis j3^] (of a case of Italian or Californian 
wine). The predictions are ro ~ bh^Mpc (Davis) and r > 50h~ 1 Mpc (for 
the volume limited samples with M < —19.5), so the threshold for the bet was 
set (by the referee N. Turok) to be the geometric average between the values: 
r = lbh^Mpc. 
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